function [S_min,fval] = PSO(T_range,lb,ub)
% PSO 粒子群算法 只能加边界约束lb,ub
fun = @(S)AvgGdop(T_range,S);
[S_min,fval] = particleswarm(fun,6,lb,ub);